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ABSTRACT 



We present an updated determination of the z ~ 4 QSO luminosity function 
(QLF), imp roving the qual i ty of the determination of the faint end of the QLF 
presented in lGlikman et al.l (120101 ). We have observed an additional 43 candidates 
from our survey sample, yielding one additional QSO at z = 4.23 and increasing 
the completeness of our spectroscopic follow-up to 48% for candidates brighter 
than R = 24 over our survey area of 3.76 deg 2 . We study the effect of using K- 
corrections to compute the rest-frame absolute magnitude at 1450A compared 
with measuring M1450 directly from the object spectra. We find a luminosity- 
dependent bias: template-based .^-corrections overestimate the luminosity of 
low-luminosity QSOs, likely due to their reliance on templates derived from higher 
luminosity QSOs. Combining our sample with bright quasars from the Sloan 
Digital Sky Survey and using spectrum-based Mi 450 for all the quasars, we fit a 
double-power-law to the binned QLF. Our best fit has a bright-end slope, a = 
3.3±0.2, and faint-end slope, j3 = 1.6^'^. Our new data revise the faint-end slope 
of the QLF down to flatter values similar to those measured at z ~ 3. The break 
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luminosity, though poorly constrained, is at M* = — 24.1^ 9 , approximately 
1 — 1.5 mag fainter than at z ~ 3. This QLF implies that QSOs account for 
about half the radiation needed to ionize the IGM at these redshifts. 



Introduction 



The primary observable that traces the evolution of QSO populations is the QSO lu- 
minosity function (QLF) as a function of redshift. The QLF can be well represented by 
a broken power-law: = <fi*[(L/L*) a + (L/ L*)°]~ l , where L* is the break luminosity 

Current measurements poorly constrain the corresponding break absolute magnitude to be 
M* 450 ~ —25 to —26 mag, at A res t = 1450A. T he bright-end slope, «, appears to evolve and 
flatten toward high redshifts, beyond z ~ 2.5 ([Richards et al.l 120061 ) . The faint-end slope, (3, 
has been measured to be aroun d —1.7 at z < 2.1; it is poorly constr ained at higher redshifts, 
but appears to flatten at z ~ 3 (IHunt et al.ll2004uSiana et al.ll2008l ). At yet higher redshifts, 
the situati on is much less clear due to the relat ively shallow flux limits of most surveys to 
date 



e.g. 



Wolf et al.ll2003l : iRichards et alJl2006h . The shape of the QLF at z > 3 is still not 



well measured, and the evolution of L* and the faint-end slope remain poorly constrained. 

We recently measured the faint end of the QLF at z ~ 4 and found an unexpected 
result: the QLF appeared to rise steeply towards faint luminosi ties, steeper than at lower 
redshifts, and in apparent conflict with some theoretical models (jGlikman et al.ll2010l . here- 
after, Paper I). This implied a more complex co-evolution of AGN and galaxies in the first 
2 Gyr of the Universe, and suggested that AGN are a more significant contributor to the 
metagal actic ioni zing UV flux than previously estimated. However, the faint-end slope com- 
puted in lPaper II relie d on the fain t est bin , which was highly incomplete with only two QSOs 
confirmed. Recently, llkeda et al.l (120101 ) measured the faint end of the QLF from z ~ 4 
QSOs in the COSMOS field and found a lower space density by a factor of ~ 3 to 4 and a 
flatter slope. The COSMOS survey covers 43% smaller area than our survey and has three 
times fewer QSOs (8 versus 24). In this paper we present new spectroscopic observations 
from our survey which increase our completeness and yield an improved estimate of the 
z ~ 4 QLF down to M1450 = —21 mag, primarily due to our improved determination of the 
contamination rates in the lowest luminosity bin used to compute the QLF. 



We use standard cosmological parameters throughout the paper: H 
Mpc" 1 , Q M = 0.30, and Q A = 0.70. 



70 km s 
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2. Sample Selection and Previous Results 



We directly measured the z = 4 QLF by taking advantage of two publicly available, 
deep, m ultiwavelength imaging surveys. The selection methodology is described in detail in 
Paper J. Our QSO s ample is derived from the NOAO Deep Wide-Field Surv ey (NDWFS; 
Jannuzi fc Deylll999[ ) and the Deep Lens Survey (DLS; IWittman et al.l 120021 ). which reach 
magnitude limits of R > 24 and provide multiwavelength imaging in B, V, R, I, and z, 
depending on the survey. We selected subfields from these surveys amounting to a total 
area of 3.76 deg 2 (1.71 deg 2 in ND WFS and 2.05 deg 2 in DL S). We created source catalogs 
from these images using SExtractor (IBertin fc Arnoutslll996l ) in dual-image mode, requiring 
a detection in the i?-band. 

To derive and understand our selection function, we built a library of model QSO 
spectra spanning the redshift range 3.5 < z < 5.2, including Monte Carlo realizations of 
Lya forest absorbers along the line-of-sight, continuum slope, and Lya equivalent width, 
to determine the location of these QSOs in color-color space. We applied color cut s to our 
catalogs, resulting in 148 QSO candidates, 30 of which are brighter than R = 23. In lPaper I 
we reported 28 spectra, concentrating on the brightest sources, of which 23 were QSOs with 
3.74 < z < 5.06. Twenty of these are brighter than R = 23. Interestingly, two QSOs at z ~ 4 
had anomalous N IV] 1486A emission, suggesting the presence of very massive, metal-poor 
star formation in their hosts (i.e., Population III; iGlikman et al.l 120071 ). 



We used the modeled QSO spectra to compute the selection probability of a QSO in our 
survey as a function of R magnitude and redshift. This allowed us to compute the volume 
density of quasars as a function of luminosity, $(M 145 o,2 = 4), using the l/V a method 
f lPage fc Carrerall200o| ). 



We computed the QLF at a fiducial wavelength that, for high-redshift, optically-selected 
QSOs, is 1450A. This wavelength lies in a part of the rest-frame ultraviolet that is free from 
emission lines, and is a proxy for the conti nuum luminosity. Lower red shift s urveys typically 
derive the QLF at rest-frame 5-band (e.g.. |Peilll995uCroom et al.ll2004l ) while iRichards et al. 
( 120061 ) use the absolute magnitude of a z = 2 quasar in the i-band, Mi(z = 2), for the SDSS 
QLF. 

A standard way to convert an observed magnitude to a monochromatic rest-frame ab- 
solute magnitude at 1450 A is to add an appropriate i^-correction, which is often computed 
eith er from a QSO template or simulated QSOs, such as our library of model QSO spectra 



sec 



Paper J ). At low redshifts, optical spectra do not cover rest-frame 1450 A and K- 



corrections are the only option. However, if rest-frame 1450 A is covered by the spectrum 
of a quasar, a more direct way to compute M 1450 is to convolve a spectrum that has been 



-4- 



shifted to the rest frame with a top hat filter centered at 1450 A, and apply a distance 
modulus. To account for slit losses, we convolve our spectra with a photometric filt er (I or 
z, depending on the survey) and scale the fluxes to the image-based photometry. In iPaper I 
we measured M 1450 for our QSOs using both methods: directly from the spectra and by ap- 
plying a i^-correction determined from our simulated model spectra. We found that M1450 
measured directly from the spectra were ~ 0.3 magnitudes fainter, on average, than Mi 450 
computed with i^-corrections. 



In IPaper 11 our spectroscopic completeness for R < 23 was 73% and we computed the 
QLF using only these brightest QSOs. We found that regardless of the method for computing 
^1450; the faint-end slope of the QLF was steeper than at lower redshifts when fit by a 
single power-law, $(M) = $*iO" a4(/3 " 1)A/ . The best-fit slopes were (3 = -1.6 ± 0.2 and 
(5 = —1.2 ± 0.2 for the .fC-correction-based and spectrum-based M 1450 , respectively. Once 
we added the three QSOs with R > 23 (whose spectroscopic completeness was only 5%) the 
QLF became even steeper with (3 = — 2.0±0.2 and (3 = — 2.5±0.2, for the .fT-correction-based 
and spectrum-based M1450, respectively. 

These steep slop es persisted when we com bined the faint QSOs w ith the QLF derived 

, and a double power- 
tood ) and the faint- 
er = -2.4 ± 0.2, 



from SDSS quasars in lRichards et al.l (120061 ) and lFontanot et ajj (12 007 



law fit to the combination of the bright-end QLF from iRichards et al. 
end points based on spectrum-based Mi 450 yielded a bright-end slope 
and faint-end slope, = —2.3 ± 0.2, without a well-constrained break luminosity; this was 
effectively a single power law. 

There were two immediate cosmo logical implications to this result: (1) models of the 
evolution of faint AGN at ~ 1 Gyr after the end of the reionization, and possibly all the 
way into the reionization era, would need to be revised; (2) AGN were a more significant 
contribut or to the metagalactic ionizing UV flux at these epochs, affecting th e evolution of 



the IGM (IHaiman et a 



200ll : IWvithe fc Loebll2003l : IShankar & Mathurl 120071 ). However, as 



emphasized in IPaper J . these results were based on only a few QSOs in the faintest bin, 
23 < R < 24. 



3. Additional Spectroscopy 

We sought to improve the measurement of the faint-end slope by obtaining additional 
spectra of R > 23 candidates. On UT 2010 May 17-18 we obtained 38 additional spectra with 
LRIS on Keck I. To maximize our efficiency, we observed using slitmasks, selecting candidates 
that lay within a single slit-mask (4' x 8'). We observed 22 candidates in NDWFS, 19 of 
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which were fainter than R = 23 magnitudes. We observed 16 DLS candidates, all fainter 
than R = 23. In addition, we obtained 5 spectra with DEIMOS on UT 2010 May 14 of 
ND WFS candidates as part of an observing program studying z ~ 4 LBGs in the Bootes 
field Lee et al.lboioh . 



Even though we increased our spectroscopic completeness by more than a factor of two, 
only one z = 4.23 QSO (R = 23.20) was found in the NDWFS field. No QSOs were found 
in the DLS field. In Figured] we show the /^-magnitude distri bution o f our candidates. On 
top we show the distribution of candidates with spectra from iPaper II in solid bins and the 
newly obtained spectra in shaded bins. The bottom panel shows the efficiency of finding 
QSOs as a function of magnitude. It is clear that our selection efficiency drops significantly 
beyond R = 23. 

The spectra that were not classified as QSOs fell into three categories: (1) featureless, 
unidentifiable spectra, (2) Lyman Break Galaxies (LBGs) at similar redshifts (z ~ 3 — 4), 
and (3) carbon stars. The objects with featureless spectra are not QSOs because they have 
smooth continua to the shortest wavelength end of the LRIS red camera (A > 5600A) and 
no lines were identified. The typical rms of our spectra is ~ 10~ 19 erg s _1 cm~ 2 A -1 , and 
the faintest Lya lines in our survey has a peak flux density of ~ 4 x 10~ 18 erg s _1 cm -2 A -1 . 
With signal-to-noise ratio of > 10, Lya has an unmistakeable signature in these spectra and 
its absence virtually guarantees that the object is not a QSO. 



LBGs at the same redshifts as the QSOs in our survey are a natural contaminant since 
their spectra also show the strong spectral break at ~ 6000A due to absorption by the Lyct 
forest. As we probe fainter fluxes we intersect the bright-end of the LBG luminosity function. 
Steidel et all (119991 ) . using a somewhat different color selection at z ~ 4, estimated a surface 
density of ~ 15 LBGs per deg 2 down to Jab — 23.5 mag and ~ 65 LBGs per deg 2 down to 
Jab — 24. This amounts to 56 and 244 LBGs expected in our area for the two magnitude 
limits, respectively, which is roughly consistent with the number of interlopers in our faintest 
bins. 



Carbon stars ar e common interlope rs in high-redshift QSO surveys f lRichards et al. 



20021 ) and vice- versa (IDownes et al.ll2004j ). A more detailed discussion of our spectroscopic 
results will be presented in Paper III. (Glikman et al., in preparation). 



4. Results 



In IPaper t we presented two versions of the QLF using absolute magnitudes computed 
from i^-corrections as well as directly from the spectra. Since the i^-corrections rely on 
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a QSO template, distribution of UV continuum slopes, and distribution of Lya equivalent 
widths, the latter two of which are modeled as Gaussian distributions, the luminosity derived 
from these i^-corrections is heavily dependent on assumptions about the spectral properties 
of QSOs. The measurement of M 1450 from spectrophotometry more directly reflects the 
true luminosity of each QSO. Therefore, going forward we only use M1450 measured from 
spectroscopy. 



We extend our QLF to higher luminosities using QSOs from SPSS. In IPaper j we 
compared our f aint-e nd QLF with the SDSS-based QLF from iRichards et al.l ( 120061 ) and 
Fontanot et al.l (120071 ). Although both bright-end QLFs used the same QSOs from SDSS, 
they derived volume densities that differed by more than a factor of 2 at their faintest bins . 
This was due to different m ethods of deriving the selection function. iFontanot et al.l (120071 ) 
used the QSO template from Cristiani fc Viol (199ol) and line and continuum properties based 
on the SDSS spectroscopic sample itself. IRichards et al.l (120061 ) . using a method more sim- 
ilar to ours, built a library of model QSO spectra with Gaussian distributions of line and 
continuum properties as we have done. Therefore, to construct a QLF that will consistently 



span t he full range of luminosities we use the selection function derived by IRichards et al. 



( 120061 ) except that we re-compute M 1450 for the SDSS QSOs directly from their spectra. 



We identified QSOs from IRichards et al.l ( 120061 . Table 5) with 3.8 < z < 5.2 to match the 
redshift range of our QSO selection; 399 QSOs were selected. We obtained their spectra from 
the SDSS database and measured M 1450 from each spectrum directly in the same manner 
as done for our f aint QSOs. Figure [2] compares M 1450 from ^-corrected magnitudes (using 
Equation 3 from Richards et al. 20061 ) with the spectrophotometrically derived values. We 
include our faint QSOs as squares (DLS) and triangles (NDWFS). We find that measuring 
M 1450 from .fT-corrections introduces a bias that overestimates the luminosities of QSOs. 
The effect appears to be more pronounced at fainter luminosities. This may be because the 
QSO template used to derive ^-corrections is composed of luminous QSOs. Luminous QSOs 
have more prominent blue bumps and therefore have stronger rest-frame UV emission than 
lower luminosity QSOs. Fainter objects have weaker blue bumps and may be overcorrected 
by the i^-correction technique, resulting in the bias seen in Figure [2j 

We use these new a bsolute magnitu d es to compute $(M) using the value of the selection 
function from Table 5 of IRichards et al.l (120061 ) and computing the available volume for each 
quasar, V a , with the limiting magnitude of % = 20.2. We bin the QLF in one-magnitude 
bins and find that the faint end of the SDSS QLF matches quite well with the bright end of 
our survey (Figure E]). Note that the faintest bin of the SDSS QLF is an outlier, li kely due 
to incompleteness. The open red circles show the QSO volume density derived in IPaper I . 
Table [1] lists the space densities of the binned QLF from SDSS and our QSOs. We compute 
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uncertainties for the bins with small numbers of QSOs (< 50) using iGehreld (119861 ) . Our 
improved spectroscopic completeness lowers the space density of the lowest luminosity bin by 
a factor of 12, or approximately 1.2 sigma relative to the poorly constrained value presented 



in 



Paper I . 



5. The Shape of the QLF 



We fit the QLF using the STY maximum likelihood method (jEfstathiou et al.l 119 88'). 
Initially we fit a single power-law function, which was the best-fit function found in lPaper I . 
The likelihood is maximized for a single power-law with a = —3.3 ± 0.1 at 68% confidence, 
±0.2 at 90% confidence, and ±0.3 at 99% confidenc e. This value i s effect ively a measure of the 
bright-end slope of the QLF and is consistent with ICroom et al.l ( 12004 ) for the 0.4 < z < 2.1 



range. 

Similarly, we measured the slope of a single power-law at the faint end, using only 
quasars fainter than M1450 = —24. The likelihood is maximized with /3 = —1.6 ±0.2 at 68% 
confidence, ±0.3, —0.4 at 90% confidence, and ±0.4, —0.6 at 99% confidence. 



We also fit the binned QLF to a double power-law function, 



$(M) 



lQ0A(a+l)(M-M*) _|_ lQ0A(/3+l)(M-M*) ' 



(1; 



omitting the outlying SDSS point at Mi 450 = —24.5. The resultant curve is plotted in 
Figure [3] with a solid black line. The shaded region represents the la uncertainties, which 
we determine by simulating 10,000 random QLFs based on our binned <£>(M) measurements 
and their associated errors. The best-fit parameters are: $* = $(M*) 



i-3+6;l x 10- 



Mpc 3 mag 



a 



-3.3 ±0.2, p = — 1.6±g;|, and M* = -24.ll?;£ (x 2 = 1.12). Figure 1 



shows the distribution of (3 vs. M* from our simulations, with the best-fit value overplotted 
(filled circle). 



The results from llkeda et al.l (120101 ) are overplotted in Figure [3] as blue squares. The 
dot-dash line is their best-fit double power-law. While our newly measured faint-end slope 
agrees with their value of (3 = —1.67+0 17) our normalization (<&*) is higher by a factor of 4 
and our break luminosity is fainter by 0.3 mag. 



We overplot in Figure |3] the QLF at z ~ 3 from ISiana et al.l (120081). who selected QSO s 
using infrared and optical photometry in the SWIRE Legacy Survey (ILonsdale et al.ll2003l ). 
The dashed line is their best-fit QLF to only the SWIRE data, while the dotted lines are fits 
to their SDSS±SWIRE QLF. Both the bright- and faint-end slopes are comparable at z ~ 3 
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and z ~ 4. The most striking evolution occurs at M* which, while poorly constrained, is 
1 — 1.5 magnitudes fainter at z ~ 4 than at z ~ 3, according to our fit. 

This is apparently at odds with the so-called "cosmic downsizing" where the space 
density of AGN peaks at different times for different luminosities. Low luminosity AGN 
are more abundant at low redshifts and the space density of high lu minosity AGN pea ks at 



earlier times. This is seen in the evol ution of the X-ray QLF (e.g., iBarger et al.ll2005l ) and 



in the optical at lower redshifts (e.g., ICroom et al.l 120091 ) . However, this has not been well 



explored past the peak of the quasar epoch 



two magnitudes deeper than iRichards et al 



2.5 -3). Based on QS Os in SDSS Stripe82 



(120061 ) , IJiang et al.l (120061 ) did not see evidence 



for AGN downsizing at these redshifts. 



6. The Contribution of Quasars to the UV Radiation Field at z ~ 4 



Following the same calculation as in iPaper J . but using the new parameters of the 
double-power-law QLF, we calculate the emissivity of quasars at 1450A to be 6 1450 = 
3.5l2 4 4 X x 10 25 erg s^ 1 Hz -1 Mpc~ 3 , ~ 50% of our earlier result. ISiana et al.l (120081 ) com- 

3.2 to be ei45o = 5.1 x 10 24 erg s" 1 



pute the specific luminosity density from QSOs at z 
Hz -1 Mpc -3 . This is a factor of ~ 7 lower than what we compute at z 



4. 



Following the formalism of iMadau et al.l (119991 ) . N\qm = 2.4 x 10 51 s _1 ionizing pho- 
tons are needed to ionize the IGM at z = 4.15 (the median redshift of our survey). Our 
parameterized QLF produc e s Aqs q = l-5it\[ x 10 51 s -1 , or ~ 60 ± 40% of the photons ion- 



izing the IGM. ISiana et al.l (120081 ) argue that at z ~ 3 QSOs contribute just under half the 
needed photon s to io nize the IGM, which is roughly consistent with our result. In addition, 
Vanzella et al.l (120101) find that z ~ 4 L BGs in the the Great Observatories Origins Deep 
Surveys (GOODS; iGiavalisco et al. 2004 ) account for < 20% of the photons needed to ionize 
the IGM, leaving QSOs as the likely dominant ionizing source at this redshift. 
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Table 1. Luminosity Function 



M 1450 Bin Center (M 1450 ) a $ ^Vqso 

(mag) (mag) (10 -8 Mpc~ 3 mag -1 ) 



SDSS 



-28.5 
-27.5 


-28.45 
-27.33 


o.2oi8:83 


2 

41 


-26.5 


-26.46 


0.93 ±0.07 


169 


-25.5 


-25.70 


4.3 ±0.5 


102 


-24.5 


-24.72 


0.418:1 


4 




NDWFS- 


-DLS 




-25.5 


-25.37 


24^ 3 


7 


-24.5 


-24.71 


o 0+8.5 
0-0-4.8 


3 


-23.5 
-22.5 


-23.47 
-22.61 


143+77 
3071^3 


7 
5 


-21.5 


-21.61 


^°^-280 


2 



a Mean magnitude of the QSOs in each bin. 
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Fig. 1. — We have improved the completeness of our spectroscopic follow-up of our QSO can- 
didates since Paper I. Top - -R -magnitu de distribution of our quasar candidates. Overplotted 
are objects with spectra from lPaper II (solid) and this work (shaded, additional objects ob- 
served May 2010). Bottom - Efficiency of the spectroscopic sample as a function of i?-band 
magnitude. Note that our efficiency drops significantly for R > 23. 



-28 -27 -26 -25 -24 

M from Spectrum 

1450 r 

Fig. 2. — The fT-corrected absolute magnitudes at rest-frame 1450A versus the M1450 mea- 
sured directly from the QSO spectra. The SDSS QSOs are represented with plus signs, 
squares are DLS QSOs, and triangles are QSOs from NDWFS. The thick line is the best fit 
relationship to all the points, with a slope of 0.9 and an offset of —3.3. The thin line shows a 
one-to-one relation. M1450 measured from spectrophotometry are significantly fainter (par- 
ticularly at lower luminosities) than those derived from ^-corrections. This may be because 
the QSO templates used were derived from high-luminosity objects, yet are applied to the 
full range of source luminosities in our sample. 
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Fig. 3. — The QSO luminosity function at z ~ 4, combining the QLF from SDSS [triangles) 
and our new, updated data {filled red circles). The points are plotted at the mean luminosity 
of each bin. We overplot the QLF from Paper I with open red circles, to show the change 
with increa sed spectroscopic c ompleteness. The blue squares are the space densities of z ~ 4 
QSOs from llkeda et all ( 120101 ) and the dot-dash line is their best-fit double power-law. The 
lower right-hand legend lists the best-fit parameters to a double power-law {solid line) shaded 
region represents the la uncertainties. Dashed and dotted lines show the z ~ 3 QLF from 



Siana et all ( 120081 ). representing the different fits to their QLF. W hile the faint-end s lopes 
agree with our measurements, our normalization, $*, is higher than llkeda et all ( 120101 ) by a 
factor of ~ 4. 
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Fig. 4. — We show the probability density of the faint-end slope, /3, vs. break absolute 
magnitude, M*, based on 10,000 Monte Carlo simulations of our binned QLF. Contour 
levels are 10%, 50%, 68% and 90% of the peak density. The best-fit parameters to the real 
QLF is plotted with a filled circle. 



